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During  the  period  covered  by  this  semi-annual  report  (1  January  to  30  June 
1975),  seven  general  topics  of  research  have  been  investigated.  I.  A st-udy  of 
anomi-lous  events,  earthquakes  which  plot  near  the  explosion  population  in  a Mg-ra^ 
diagram,  has  concluded  that  all  such  anomalous  events  in  the  eastern  Hinalayas  are 
natural  earthquakes.  II.  A study  of  a foreshock-mainshock -aftershock  sequence  is 
being  conducted  to  determine  if  there  are  factors  other  than  source  depth  focal 
mechanism  that  control  the  surface  wave  magnitude  of  small  events.  The  study  indi- 
cates that  finiteness  of  the  source  can  play  a role  in  controlling  for  events 
smaller  than  m^^  5.2.  III.  Three  studies  are  in  progress  on  the  topic  of  Intra- 
plate Tectonics:  a)  focal  mechanisms  for  Intraplate  earthquakes  in  Central  Asia 
by  IWNSJ  and  mainland  Chinese  data;  b)  analysis  of  an  accelerogram  of  an  intra- 
plate  earthquake.  Blue  Mtn.  Lake,  N.Y.;  md  c)  In  Situ  strain  relief  measurements  in 
an  intraplate  area,  the  Adirondack  Mts.,  'l.Y.  IV.  A study  comparing  the  motions  in 
a half  space  near  a fault  v?lth  those  which  occur  in  a whole  space  if  the  surface 
were  abse.  t has  been  completed.  The  results  are  t'-at  the  motion  on  the  half  space 
may  be  approximated  by  doubling  those  in  an  infinite  space  when  the  angle  of  inci- 
dence is  less  than  30°.  V.  A study  is  currently  in  progress  t(  test  the  validity 
of  various  mouels  of  the  core-mantle  boundary  using  seismograms  synthesized  by  the 
Langer  method  of  frequency  domain  calculations.  VI.  In  another  study,  the  excita- 
tion of  higher  mode  surface  waves  by  shallow  events  is  being  exairdned  . VII.  A 
study  is  also  under  way  on  the  possible  interrelation  among  bouyant  zones,  great 
earthquakes,  and  dynamic  boundaries  of  subductlon. 
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3. 


MAJOR  SCIENTIFIC  ACCOMPLISHMENTS 


In  the  following  paragraphs,  scientific  accomplishments  for  January  through 
June  1975  are  summarized  following  the  Itemization  In  the  statement  of  work  of 
this  contract. 


Line  - Item  0001b 


1.  Discrimination  of  Anomalous  Events  In  the  Eastern  Himalayas 

Ongoing  research  toward  the  discrimination  of  anomalous  events — those 
natural  earthquakes  that  plot  in,  or  near,  the  explosion  population  on  an 
Mg-m|^  diagram — has  led  to  the  Identification  of  all  such  events  located  near 
the  eastern  Himalayas  as  natural  earthquakes.  Most  of  the  reported  anomalous 
events  (about  25)  are  located  In  Tibet  in  a particular  tectonic  setting 
associated  with  the  eastern  termination  of  the  Himalayas.  At  this  point 
several  tectonic  trends  of  the  Himalayas,  an  area  of  continent-continent 
collision,  converge  toward  an  Intersection  with  the  northward  projection  of 
the  Andaman-Burman  arc.  Thus,  the  tectonic  regime  of  these  anomalous  events 
Is  not  part  of  the  Himalayas  but  rather  appears  to  be  related  to  "end-effects" 
at  their  eastern  terminus.  Previous  investigators  have  suggested  that  anomalous 
events  result  from  strong  excitation  and  efficient  propagation  of  body  waves. 

Our  study  of  the  regional  tectonics  of  the  area  reveals  several  potential 
sources  of  poor  surface  wave  excitation.  The  crust  of  the  Tibetan  Plateau, 
which  from  meager  data  appears  to  be  almost  twice  as  thick  as  normal  continental 
crust,  offers  a medium  for  the  occurrence  of  earthquakes  at  moderate  focal 
depths.  Rayleigh  waves  of  the  fundamental  mode  are  poorly  excited  by  earthquakes 
at  such  ('epths.  The  thick  crust  also  provides  a wave  guide  for  the  efficient 
propagation  of  Love  waves  and  higher-mode  Rayleigh  waves,  and  it  Is  found 
that  observations  of  these  waves  provide  a discriminant  for  some  of  the  anomalous 
events.  In  particular  observations  of  higher  modes  offer  an  especially 
attractive  dlscriml'.ant . Since  the  degree  of  excitation  of  various  modes 
depends  mainly  on  focal  depth,  no  explosion  data  are  required  for  comparison 
of  earthquakes  and  explosions.  These  higher-mode  observations,  combined 
with  more  reliable  m^^  determinations,  Identify  all  of  the  reported  anomalous 
events  In  Tibet — with  the  exception  of  a single  sequence  of  events — as  natnral 
earthquakes.  The  remaining  sequence  of  anomalous  events,  none  of  which  is  in 
the  explosion  population,  occurred  xn  a very  limited  area  (about  50  km  or  less 
on  a side)  where  a concentration  of  tectonic  stress  may  be  anticipated.  Thus, 
these  few  remaining  anomalous  earthquakes  may  be  related  to  high  tectonic 
stress,  perhaps  In  association  with  the  formation  of  a new  fault.  Hlgh-stress 
earthquakes  often  have  small  source  dimensions,  and  hence  only  weak  surface 
waves  are  excited.  This  study  places  the  occurrence  of  these  remaining 
anomalous  events  with  the  complex  tectonic  setting  of  the  eastern  terminus 
of  the  Himalayas.  Observations  of  dllatatlonal  first  motions,  however,  allow 
these  events  .:o  be  discriminated.  Thus,  all  of  the  reported  anomalous  events 
In  Tibet  are  Identified  as  natural  earthquakes.  To  our  knowledge,  all  events 
from  Tibet  or  elsewhere  that  have  been  discussed  In  the  literature  as  being 
anomalous  can.  In  fact,  be  Identified  as  earthquakes  using  one  or  more  of  the 
above  mentioned  techniques. 
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Higher  Mode  Excitation  by  Phallow  Earthquakes 

The  excitation  of  higher  inodes  by  shallow  events  large  enough  to  allow 
focal  mechanism  determination  is  being  examined.  The  spectra  of  the  funda- 
mental mode  Love  and  Rayleigh  waves  observed  at  several  stations  are  used  to 
establish  the  depth  of  the  events.  Since  all  of  the  events  studied  are  in  a 
limited  geographical  area,  a comparative  event  technique  can  be  used  to  remove 
propagation  path  effects  from  the  spectra.  Once  the  depths  and  mechanisms 
are  known,  the  events  can  be  used  as  standards  for  establishing  the  detailed 
characteristics  of  excitation  and  propagation  of  the  higher  modes.  In  turn, 
these  characteristics  can  be  helpful  In  classifying  other  seismic  sources, 
particularly  those  small  events  for  which  higher  modes  are  observed,  but  which 
are  too  small  for  reliable  measurement  of  the  complete  spectra  of  the  funda- 
mental mode.  For  some  small,  anomalous  events  In  Kazakh,  an  M based  on  the 
higher  mode  amplitudes  may  be  more  reliable  than  an  M based  on  the  short- 
period  fundamental  mode  Rayleigh  waves.  ® 


Line  - Item  OOOld 


Causes  of  Variations  Within  a Centra 1-Asian  Earthquake  Sequence 


A foreshock-malnshock-af ter shock  sequence  which  recurred  in  June,  1971 
In  the  Kirglz-Kinkiarg  border  region  is  being  studied  to  determine  what  factors 
other  than  source  depth  and  focal  mechanism  control  the  surface  wave  magnitude 
of  small  events.  The  National  Earthquake  Information  Center  (NEIC)  assigned 
a body  wave  magnitude  of  5.6  to  both  the  foreshock  and  malnshock,  but  the 
surface  waves  generated  by  the  malnshock  were  approximately  5 times  longer 
than  the  foreshock.  ISC  assigned  m^^  values  5.4  and  5.3  to  the  foreshock  and 
main  sheck,  respectively.  The  surface  wave  magnitude  could  be  affected  by 
differences  in  depth,  mechanism,  stress  drop  or  finiteness,  but  all  possible 
path  and  instrument  phenomena  are  eliminated  by  choosing  a sequence  which  Is 
limited  in  duration  and  location.  The  mechanism  and  depth  of  the  events  can 
be  deduced  from  comparison  of  the  azimuthal  variations  in  the  frequency- 
dependent  radiation  patterns  of  the  amplitude  and  Initial  phase  of  the 
Rayleigh  and  Love  waves. 

Evernden  (1975)  suggested  that  the  primary  factors  controlling  M of  small 
events  (less  than  about  mj,  5.5)  are  depth  and  focal  mechanisms  and  t^at 
stress  drop  and  source  length  are  unimportant.  However,  differences  in 
mechanism  and  depth  do  not  appear  to  adequately  explain  the  dramatic  varia- 
tion in  surface  wave  amplitM/»c  within  the  series  of  aftershocks.  For  example, 
one  aftershock,  with  m.  5.2,  generated  surface  waves  roughly  12  times  larger 
than  another  aftershock  with  mj,  5.1.  The  contrast  cannot  be  attributed  to 
an  unusual  combination  of  depth  and  mechanism  effects  on  the  surface  wave 
amplitudes,  because  the  long-period  body  waves  are  affected  to  a similar 
degree.  The  earthquakes  in  the  sequence  can  be  divided  into  two  categories* 
those  which  excite  long-period  signals  efficiently,  such  as  the  main  shock  and 
several  aftershocks;  and  those  which  do  net,  such  as  the  foreshock  and  the 
majority  of  aftershocks.  The  primary  difference  between  the  two  groups  may 
be  In  the  finite  size  of  the  source.  As  observed  teleseismlcally  the  second 
group  typically  generates  a very  simple,  short-period  f’-wave  with  most  of  the 
energy  concentrated  within  the  first  1 to  2 sec.,  suggesting  nearly  a point 
source.  However,  earthquakes  within  the  first  group  generate  a much  more 
complex  P-wave,  with  substantial  energy  spread  out  o^’er  the  first  10  to  15 
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sec.  These  observations  suggest  that  flnltaness  of  the  source  can  play  an 
important  role  in  controlling  Mg  : even  for  earthquakes  with  body  wave 

magnitude  5.2  or  smaller. 


Line  - Item  OOOle 

1.  Analysis  of  an  Accelerogram  of  an  Tntraplate  Earthquake 

In  advancing  our  understanding  of  possible  driving  mechanisms  of  plate 
tectonics,  Sykes  & Sbar  (1973,  1974)  have  considered  the  focal  mechanism  of 
intraplate  earthquakes.  Most  Intraplate  earthquakes  were  interpreted  to 
be  of  thrusting  mechanisms,  suggesting  horizontal  compression  as  the  dominant 
component  of  stress.  In  fact,  Sbar  and  Sykes  (1973,  1174)  have  suggested 
much  of  eastern  North  America  Is  in  a state  of  high  horizontal  compressive 
stress.  We  have  recently  completed  a detailed  study  of  one  Intraplate 
earthquake  In  the  Eastern  United  States.  The  principal  results  are  that  an 
accelerogram  obtained  at  Blue  Mt.  Lake  is  remarkable  for  the  simplicity  of 
its  S-wave  pulse.  This  is  due  to  1)  a nearly  complete  absence  of  scattering 
and  reflections  as  second  arrivals  on  the  accelerogram  and  2)  a very  elemen- 
tary earthquake  source.  The  earthquake  identified  with  this  accelerogram 
had  a magnitude  m^j  = 2.2  and  a hypocentral  distance  of  about  1 km  from  the 
accelerometer.  Analysis  of  the  S-wave  indicates  the  earthquake  had  a moment 
of  5 to  6.3  X 10^®  dyne-cm,  and  a source  radius  of  20  to  40  m.  When  the 
accelerogram  is  integrated  to  obtain  displacement,  there  is  a step  offset  of 
about  5 microns  associated  with  a near  field  component  of  the  S-wave  pulse. 

The  S-wave,  including  the  step  offset,  can  be  matched  in  remarkable  detail 
by  a dislocation  model. 

2.  Near  Field  Displacements 

Vie  have  made  a systematic  comparison  of  motions  on  a half  space  near 
a fault  with  those  which  would  occur  in  a whole  space  if  the  surface  were 
absent.  We  use  the  Green's  function  for  a half  space  given  by  Johnson  (1974). 
Use  of  this  Green's  function  allows  four  major  surface  effects  to  be  studied: 
the  amplification  of  all  waves;  the  phase  shift  of  SV  waves  incident  at 
angles  more  grazing  than  critical;  the  SP-phase  which  is  a P-wave  converted 
from  an  SV  wave  incident  at  the  critical  angle;  and  Rayleigh  waves.  Our 
objective  is  a critical  evaluation  of  the  assumption,  used  frequently  in 
dislocation  modeling,  that  the  free  surface  can  be  reasonably  accounted  for 
by  doubling  the  whole  space  amplitude  (eg.  Kanamori^  1972;  Trlfunac,  1974; 
Trlfunac  and  Udwadla,  1973;  Anderson,  1974). 

In  doing  this,  we  have  developed  a practical  way  of  computing  displace- 
ments on  the  surface  of  a half  space  as  a two  step  process.  The  first  step 
is  calculating  and  storing  the  Green's  function  for  a given  station  for 
several  points  on  the  fault;  and  the  second  step  is  convolving  with  the 
source  time  function.  This  method  can  be  faster  than  the  corresponding 
whole  space  calculations  (which  recompute  the  Green's  function)  when  several 
models  are  tried  for  the  same  geometry,  as  in  a studying  a particular  ac- 
celerogram record.  This  method  also  gives  a clear  idea  of  what  frequencies 
are  significant  in  a given  record. 


6. 


The  results  are  that  motions  on  the  surface  of  a half  space  may  be 
approximated  by  doubling  the  amplitude  of  motion  in  an  infinite  space  when 
angles  of  incidence  are  less  than  30  degrees.  For  greater  angles  of  incidence, 
the  faulting  parameters  derived  from  a whole  space  model  in  this  manner  are 
potentially  somewhat  mislead Ing . The  P-  and  SV-  waves  from  a point  source 
in  the  half  space  can  generally  be  obtained  from  the  whole  sp<i>-a  counter- 
parts by  applying  the  appropriate  plane  wave  surface  correction.  The  SV- 
waves  for  angles  somewhat  greater  than  critical  cannot  be  derived  in  this 
way,  neither  can  the  AP-phase  of  the  Rayleigh  wive. 


Line  - Item  0001 f 


Earth  Modeling  using  Synthetic  Seismograms 

A study  is  currently  in  progress  to  test  the  validity  of  various  models 
of  the  core-mantle  boundary  using  synthetic  seismograms.  The  Langer  method 
of  frequency  domain  calculations  (Richards,  1975)  has  distinct  advantages  over 
other  methods  in  describing  the  effects  of  earth  structure  on  waves  propagating 
through  the  earth.  The  method  is  uniformly  asymptotic  in  frequency,  as  opposed 
to  virtually  all  other  approximating  methods  which  are  non-uniform  and  fall 
at  turning  points.  In  addition,  for  a ray  with  a turning  point  near  a dis- 
continuity in  the  earth,  the  Langer  approach  corrects  the  reflection  and 
transmission  coefficients  for  dependence  on  frequency,  radius  or  curvature  of 
a discontinuity  and  the  effect  of  earth  structure  near  the  discontinuity. 
Indeed,  the  observation  of  multiple  core  phases  at  distances  beyond  the  range 
predicted  by  classical  ray  theory  for  currently  acceptable  earth  models  can  now 
be  explained  in  terms  of  tunnelling  phenomenon  (Richards,  1973,  1975)  in 
which  rays  with  turning  points  .lust  above  the  core-mantle  boundary  can  leak 
energy  into  the  earth's  core. 

The  Langer  approach  is  currently  being  aoplled  to  the  synthesis  of  the- 
oretical seismograms  of  SKS  and  5KKS  and  recorded  on  long  period  instruments. 
The  method  also  easily  accounts  for  diffracted  energy  associated  with  SKS 
in  a body  wave  called  SP^^j^ffKS  (King  and  Muller,  1975).  The  behavior  of  the 
long  period  content  of  these  body  waves  should  deviate  from  classical  ray 
theory.  The  extent  and  sensitivity  of  this  deviation  for  various  earth  models 
will  be  Investigated.  The  comparison  of  synthesized  seismograms  with  data 
from  WWSSN  and  other  seismic  stations  will  serve  as  a guide  to  the  validity 
of  various  earth  models  at  the  core-mantle  transition  zone. 


Line  - Item  OOOlg 


A Study  of  Focal  Mechanisms  for  Intraplate  Earthquakes  in  Central  Asia 

A focal  mechanism  study  of  the  contemporary  tectonics  of  central  Asia  is 
being  undertaken  using  data  trom  hoth  WWSSN  film  chi])s  and  Information  from 
seismograms  recorded  in  mainland  China.  The  cpic'.ntral  distribution  in  this 
area  shows  that  there  is  no  sharp  boundary  between  the  Indian  and  Eurasian 
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plates  as  would  be  expected  If  the  contemporary  tectonics  were  simply  the 
result  of  the  intersection  of  the  two  ates.  Molnar  et  al.  (1973)  proposed 
that  stable  blocks  within  the  Eurasian  plate  may  play  an  important  role  in 
the  explanation  of  the  tectonics.  There  presented  no  quantitative  evidence 
for  this  proposal. 

The  quantitative  evidence  necessary  to  study  this  tectonic  problem  has 
not  been  available  until  the  present.  Data  from  the  close-in  Chinese  stations 
are  necessary  to  determine  reliable  focal  mechanisms  for  the  thrust  fault 
type  earthquakes  which  predominate  In  the  fold  belts  of  central  Asia.  The 
WWSSN  station  distribution  does  not  lend  itself  to  reliable  focal  mechanisms 
for  this  area. 

Through  the  analysis  of  the  combined  WWSSN  and  Chinese  data  sets,  we 
hopr;  to  obtain  more  accurate  orientations  for  the  slip  vectors  and  the 
directions  of  the  three  principal  stress  axes.  The  former  can  be  used  to 
define  the  direction  of  relative  motion  between  two  of  the  plates  and  in 
effect  aid  in  the  determination  of  the  plate  boundaries.  From  the  stress 
orientations,  we  hope  to  obtain  an  improved  understanding  of  the  relations 
between  the  seismicity  and  the  tectonic  features. 

In  Situ  Strain  Relief  Measurements  in  an  Intraplate  Area 

Twelve  strain  relief  measurements  were  made  at  5 sites  on  Potsdam  Sand- 
stone northeast  of  the  Adirondactc  Mountains,  New  York  to  test  Sbar  and 
Sykes  (1973)  suggestion  that  the  upper  crust  in  a large  area  of  eastern  North 
America  is  currently  in  a state  of  east-west  compression.  Iii  situ  strain 
was  relieved  by  overcoring,  and  detected  using  strain  gages  bonded  to  the 
sandstone  at  the  surface.  The  average  of  nine  measurements  at  four  sites  in 
the  Late  Cambrian  Keesevlllt  Member  of  the  Potsdam  Sandstone  indicate  that 
the  direction  of  maximum  expansion  upon  overcoring  is  N78°W.  These  nine 
measurements  also  indicate  that  stress  can  be  consistent  in  orientation 
over  areas  of  tens  of  square  kilometers.  The  average  of  three  measurements 
at  one  site  in  the  Nicholvllle  Member  11  km  from  the  other  four  sites 
indicate  a maximum  expansion  oriented  N18°E.  Double  overbore  tests  suggest 
that  a sizeable  component  of  the  measured  strain  is  residual.  Wv  suggest 
that  the  residual  strain  measured  in  the  Nicholvllle  Member  may  oe  related 
to  a late  Proterozoic  (Hadrynian)  or  early  Cambrian  stress  field. 

Line  - Item  QOOlh 

A Study  of  the  Relationship  among  Bouyant  Zones,  Great  Earthquakes  and 

Dynaulc  Boundaries  of  Subductlon 

The  distribution  of  large  shallow  earthquakes  along  subductlon  boundaries 
is  in  serious  disagreement  with  the  distribution  pattern  that  might  be  pre- 
dicted from  a simple  model  of  plate  tectunlcs.  That  is,  along  extensive  sections 
of  some  island  ercs  large  shocks  occurred  infrequently  or  not  at  all  during 
recorded  history.  Most  of  these  zones  of  long-term  quiescence  are  nearly 
coterminous  with  segmenl s of  the  margin  where  zones  of  seamounts,  aseismlc  ridges 
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or  other  bathymetric  highs  of  the  underthrust  slab  appear  to  be  interacting 
with  the  subduct  Ion  process.  Because  of  this  widespread  spatial  correlation 
and  because  all  or  nearly  all  elevated  portions  of  the  sea  floor  have  crustal 
roots,  aseismic  ridges  or  other  uplifted  regions  may  reveal  zones  of  relative 
buoyancy  within  oceanic  plates  which  resist  subduction  upon  collision  with  an 
active  trench.  Thus,  in  place  of  the  typical  subduction  process  at  locatlors 
where  buoyant  zones  Interact  with  the  subduction  process,  rigid  plate  motion 
may  be  preser\'ed  by  a variety  of  tectonic  responses  including:  arc-polarity 

reversal,  gradual  or  abrupt  shifts  in  the  plate  margin,  the  development  of 
cusps  between  arcs  (as  suggested  by  Vogt,  1973),  development  of  broad  zones 
of  deformation,  partial  subduction  (descent  to  pertiaps  50-100  km  but  no  farther) 
or,  possibly,  termination  of  subduction.  If  our  interpretation  is  valid,  then 
it  may  be  possible  to  explain  in  a relatively  straightforward  manner  othev 
modifications  of  the  subduction  process  such  as  gaps  in  the  line  of  active 
volcanoes,  the  presence  of  detached  slabs  of  lithosphere  in  certain  parts  of 
subduction  zones,  abrupt  changes  in  strike  and  dip  of  the  descending  litho- 
spheric plate,  the  absence  of  intermediate-depth  earthquakes  in  certain  loca- 
tions and  the  emplacement  of  ophlolite  sequences. 

In  place  of  ’’typical  oceanic"  lithosphere,  therefore,  there  may  exist  a 
broad  spectrum  of  average  densities  for  oceanic  lithosphere  and  the  relative 
buoyancies  of  the  two  slabs  near  the  subducting  margin  appears  to  be  a dominant 
influence  in  the  development  of  subduction  tectonics  and  in  the  locations  of 
great  earthquakes. 
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